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Abstract. The paper presents two algorithms for finding irreducible decomposi- 
tion of monomial ideals. The first one is recursive, derived from staircase structures 
of monomial ideals. This algorithm has a good performance for highly non-generic 
monomial ideals. The second one is an incremental algorithm, which computes 
decompositions of ideals by adding one generator at a time. Our analysis shows 
that the second algorithm is more efficient than the first one for generic mono- 
mial ideals. Furthermore, the time complexity of the second algorithm is at most 
0{n^p£) where n is the number of variables, p is the number of minimal genera- 
tors and £ is the number of irreducible components. Another novelty of the second 
algorithm is that, for generic monomial ideals, the intermediate storage is always 
bounded by the final output size which may be exponential in the input size. 



1. Introduction 

Monomial ideals provide ubiquitous links between combinatorics and commuta- 
tive algebra [24, 16]. Though simple they carry plentiful algebraic and geometric 
information of general ideals. Our interest in monomial ideals is motivated by a 
paper of [9], where they studied the connection between the structure of monomial 
basis and the geometric structure of the solution sets of zero- dimensional polyno- 
mial ideals. Irreducible decomposition of monomial ideals is a basic computational 
problem and it finds applications in several areas, ranging from pure mathematics 
to computational biology, see for example [12] for computing integer programming 
gaps, [3] for computing tropical convex hulls, [22] for finding the joins and secant 
varieties of monomial ideals, [2] for partition of a simplicial complex, [19] for solving 
the Frobenius problem, and [13] for modeling gene networks. 

We are interested in efficient algorithms for computing irreducible decomposition 
of monomial ideals. There are a variety of algorithms available in the literature. 
The so-called splitting algorithm: Algorithm 3.1.2 in [23] is not efficient on large 
scale monomial ideals. [17] gives two algorithms: one is based on Alexander duality 
[14], and the other is based on Scarf complex [4]. [18] improves the Scarf complex 
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method by a factor of up to more than 1000. Recently, [20] proposed several slicing 
algorithms based on various strategies. 

Our goals in this paper are to study the structure of monomial ideals and present 
two new algorithms for irreducible decomposition. We first observe some stair- 
case structural properties of monomial bases in Section 4. The recursive algorithm 
presented in Section 5 is based on these properties, which allow decomposition of 
monomial ideals recursively from lower to higher dimensions. This algorithm was 
presented as posters in ISSAC 2005 and in the workshop on Algorithms in Algebraic 
Geometry at IMA in 2006. Our algorithm was recently generalized by [20] where 
several cutting strategies were developed and our algorithm corresponds to the min- 
imum strategy there. Also, the computational experiments there shows that our 
algorithm has good performance for most cases, especially for highly non-generic 
monomial ideals. 

Our second algorithm is presented in Section 6. It can be viewed as an improved 
Alexander dual method ([14, 17]). It is incremental based on some distribution 
rules for "+" and "fl" operations of monomial ideals. We maintain an output list 
of irreducible components, and at each step we add one generator and update the 
output list. In [17], there is no specific criterion for selecting candidates that need to 
be updated, and the updating process is inefficient too. Our algorithm avoids these 
two deficiencies. Our analysis in Section 7 shows that the second algorithm works 
more efficiently than the first algorithm for generic monomial ideals. We prove that, 
for generic monomial ideals, the intermediate storage size (ie. number of irreducible 
components at each stage) is always bounded by the final output size, provided 
that the generators are added in lex order. This enables us to show that the time 
complexity of the second algorithm is at most 0{n'^pi) where n is the number of 
variables, p is the number of minimal generators and i is the number of irreducible 
components. 

In Section 2, we present some notations and introductory materials on monomial 
ideals. In Section 3 we discuss tree representations and operations of monomial 
ideals. 



2. Monomial Ideals 

We refer the reader to the books of [5] for background in algebraic geometry 
and commutative algebra, and to the monograph [16] for monomial ideals and their 
combinatorial properties. 

Let IK be a field and ]K[X], the polynomial ring over K in indeterminates X = 
Xi, . . . , Xn- For a vector a = (ai, . . . , a„) G N", where N = {0, 1,2,.. .} denotes the 
set of nonnegative integers, we set 



X" = 



An 

n 1 
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which is called a monomial. Thus monomials in n variables are in 1 — 1 correspon- 
dence with vectors in N". Suppose a = (ai, . . . , a^) and /5 = . . . , 6„) are two 
vectors in N", we say 

a < (3 ii Qj < bj for all 1 < j < n. 

This defines a partial order on N", which corresponds to division order for monomials 
since x°'\x^ if and only if a < We say 

a < P if a < (3 but a ^ (5. 

Also we define 

a -< (3 ii Qj < bj for all 1 < j < n. 
Then a (5 means that > bj for at least one j. 

An ideal I C K.[X] is called a monomial ideal if it is generated by monomials. 
Dickson's Lemma states that every monomial ideal in ]K[X] has a unique minimal 
set of monomial generators, and this set is finite. Denote this set to be Min(J), that 
is, 

Min(J) = {X" G / : there is no X'^ e I such that /3 < a}. 

A monomial ideal / is called Artinian if / contains a power of each variable, or 
equivalently, if the quotient ring K[X]/J has finite dimension as vector space over 
K. For convenience of notations, we define 

= 0, l<i<n. 

By adding infinity power of variables if necessary, a non- Artinian monomial ideal 
can be treated like an Artinian monomial ideal. For example, / = (x'^y^) = 
{x°°, x'^y^, y°°). Instead of adding infinity powers, we can also add powers x^" where 
Cj is a sufficiently large integer, say larger than the largest degree of Xi in all the 
monomials in Min(J). Then the irreducible components of the original ideal are in 
1-1 correspondence to those of the modified Artinian ideal; See Exercise 5.8 in [16] 
or Proposition 3 in [20]. In our algorithms belows, we will use infinity powers, but 
in the proofs of all the results, we will use powers 

An ideal J C K.[X] is called irreducible if it can not be expressed as the inter- 
section of two strictly larger ideals in ]K[X]. That is, J = Ji fl J2 implies that J = Ji 
or J = J2- A monomial ideal I is irreducible if and only if I is of the form 

for some vector (3 = . . . , 6„) G W" where N = N U {00} \ {0}. Thus irreducible 
monomial ideals are in 1-1 correspondence with /3 G N . 

An irreducible decomposition of a monomial ideal I is an expression of the 
form 

/ = n ■ ■ ■ n m^-- (1) 

where . . . , (3r ^ ■ Since the polynomial ring ]K[X] is Noetherian, every ideal 
can be written as irredundant intersection of irreducible ideals. Such an intersection 
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is not unique for a general ideal, but unique for a monomial ideal. We say that 
the irreducible decomposition (1) is irredundant if none of the components can be 
dropped from the right hand side. If (1) is irredundant, then the ideals m^^, . . . , m^*" 
are called irreducible components of /. We denote by Irr(J) the set of exponents 
of irreducible components of /, that is, 

Irr(/) = {/5i,...,/3,}. 

By this notation, we have 

1= f] mf. 

/3elrr(7) 

Note that, for two vectors a and /5, 

X" G if and only ii a j3, 

and 

m" C if and only ii f3 < a. 

A monomial ideal / is called generic if no variable Xi appears with the same 
non-zero exponent in two distinct minimal generators of /. This definition comes 
from [4]. For example, 

is generic, but 

is non-generic, as z"^ appears in two generators. Loosely speaking, we can say I2 is 
nearly generic, but 

h = {xy,yz,xz,z'^) 

is highly non-generic. Previous algorithms [17, 18] behave very different for generic 
monomial ideals and highly non-generic monomial ideals. For example, the Scarf 
complex method works more efficient when dealing with generic monomial ideals 
[17]. 

In the following sections, we always assume that we are given the minimal gen- 
erating set of a monomial ideal. Though our algorithms work for monomial ideals 
given by an arbitrary set of generators, it will be more efficient if the generators are 
made minimal first. 

3. Tree Representation and Operations 

Note that monomials are represented by vectors in N" and irreducible components 
are represented by vectors in N . To efficiently represent a collect of vectors, we use 
a tree structure. This is used in [9, 17]. This data structure is also widely used in 
computer science, where it is called a trie. 
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Tree representation. First we want to define the orderings on N" or N . Suppose 
a = (ai, . . . , a„) and (3 = (61, . . . , 6„) are two vectors in N" or N , and the variable 
ordering is Xi < ■ ■ ■ < a:„ in K[X]. We say a <iex P if aj = bj for k + 1 < j < n, 
but ttk < bk for some 1 < k < n. 

Next, suppose 5* C is a set of vectors corresponding to the generators of a 
monomial ideal / C ]K[X]. We represent as a rooted tree T of height n in a natural 
way. The tree should have IS"! leaves and the unique path of the tree from the root 
to a leaf represents a vector in S. Precisely, to represent a vector a = (ai, . . . , a„), 
we label all the nodes except the root of the path simply by a„, . . . , ai in the order 
from the root to the leaf. We regard the root as being at height 0. For two vectors 
a = (ai, . . . , ttn) and /? = (61, . . . , if aj = bj ior k + 1 < j < n but 7^ b^, then 
a and f3 share their corresponding path until height n — k. After that their children 
are listed in increasing order with respect to their coordinates. Figure 1 is the tree 
representation for I = (x^, |/^, x^y'^z'^, xy^z^, x^yz^) with variable order x < y < z. 




Figure 1. An example of tree representation. 

The tree representation for a set of irreducible components could be constructed 
in a similar manner. 

To perform the operations on sets of vectors, we need only perform on trees. We 
need three basic tree operations: Merge, MinMerge and MaxMerge. 

Merge. Given q rooted trees 7i, . . . , 7^ with the same height, merge them to form 
one rooted tree with the same height. Here we simply put the paths from all the 
trees together with repetition ignored (actually no repeated paths occur in our al- 
gorithms). We stress that no reduction work is performed under this operation. 

MinMerge. We use MinMerge(7i, . . . , Tq) to represent the set of minimal elements 
in Merge(TL, . . . , Tq). For two vectors a, [3 in Merge(TL, . . . , Tq), if a < /3, ie. x"\x^ , 
then the path for (5 should be removed in this operation. The purpose is to find the 
minimal generating set for the ideal I1 + ■ ■ ■ + Iq where % is the tree representation 
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for /j. 

MaxMerge. Similarly, the set of maximal elements in Merge(7^, . . . , 7^) is rep- 
resented by MaxMerge(7^, . . . If a < P, ie. m'^ C m", then the path for a 
should be removed in this operation. Hence, if % represents the set of irreducible 
components of Jj, 1 < z < q, then MaxMerge(TL, . . . ,7^) represents the the set of 
irreducible components of the ideal Ji fl ■ ■ ■ fl /g. 

4. Structure Properties of Monomial Bases 

In the results and their proofs below, we explicitly assume that all the ideals are 
Artinian, adding large powers if necessary where is an integer, though infinity 
powers will be used in the Algorithms and Examples. 

The monomial basis B(J) for a monomial ideal / is defined as 

B(/) = {7 G N" : ^ J}, 

which form a linear basis for the quotient ring K[X]// over K. Thus, for 7 G N", 
7 G B{I) if and only if a ^ 7 for every a G Min(J). Note that B(/) is a 6-set, that 
is, if 7 G B(J) and < 7, then n G B(J). The next lemma characterizes B[I) in 
terms of Irr(J). 

Lemma 1. For 7 G N", 7 G B{I) if and only if '~f -< jS for some (3 G Irr{I). 

Proof. Since / = ■• have X'^ G / if and only if X"' G m^, ie., 7 -/{ /?, 

for each (3 G Irr(/). Hence X'^ ^ / if and only if 7 -< /? for some (3 G Irr(/), as 
desired. □ 

We now want to express Irr(J) in terms of B(J). Since / is Artinian, for (3 = 
(61, . . . ,bn) G Irr(J), we have bi > for 1 < i < n. Define 

/5el = (6i-l,62-l,...,&n-l). 

Lemma 1 implies that, for each P G Irr(J), we have /3 1 G B(J). 
A vector 7 G N" is called maximal in B(/) if 

7 G B(J) and there is no /i G B(J) such that /i > 7- 

Lemma 2. For any vector /5 G N", /5 G Irr{I) if and only if (3 Q 1 is maximal in 
B{I). 

Proof. By Lemma 1, /?0l G B(/) if and only if there is a G Irr(/) such that /90l -< a. 
Notice that a 1 G B(/) and /3 1 -< a is equivalent to say /5 1 < a 1. Hence 
/? 1 is maximal in B{I) if and only if /5 1 = a 1, that is, /5 = a G Irr(/). □ 

The staircase diagram will help us visualize the structural properties of monomial 
ideals. For example. Figure 2 is the staircase diagram for the monomial ideal / = 
(x^, y^, x^y'^z^, xy^z^, x'^yz^). In this figure the gray points are in 1-1 correspondence 
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Figure 2. An example of staircase diagram. 

with the minimal generators, while the white points are in 1-1 correspondence with 
the irreducible components of /. Geometrically, B(/) is exactly the set of interior 
integral points of the solid. 

5. Recursive Algorithm 

For bivariate monomial ideals, irreducible decomposition is simple [15]. Suppose 

Min(/) = {x''\x''^y''\...,x''^-'y^''-\y^''} 

where ai > ■ ■ ■ > ap„i > 0, < b2 < ■ ■ ■ < bp, and Oi or bp can be infinity. Then the 
irreducible decomposition of / is 

/ = n {x^^y^') n ■ ■ ■ n n {x'^^-^y^^). 

Our recursive algorithm is a generalization of the above observation to higher 
dimensions. Let / C K[a;i, . . . ,Xn] be a monomial ideal. Suppose all the distinct 
degrees of Xn in Min(/) are 

= do < di < . . . < dg. 

For example, in / = (x'^y^) = {x°°,x'^y^,y°°), the distinct degrees in y are do = 
0, (ii = 3 and d^ = oo. We collect the coefficients of m G Min(J) as polynomials in 
Xn- Precisely, for < A; < s, let 

Ik = {CoeSx„{m) : m G Min(/) and deg^.^^ m < d^) C K[a;i, . . . , Xn-i]- 

Then 

/o £ /i £ ■ ■ ■ C (2) 

By (2), it follows that 

B{Io)DB{h)2---2B{Is). 



8 



SHUHONG GAO AND MINGFU ZHU 



For the example with I = x^y^, ?/°°), Jq = {x°°) = {0}, Ji = = (x^), 

and h = 1) = (1) = K[x]. 

We show how to read off the irreducible components of I from those of J^'s, which 
have one less variables. For any vector fi = {ui, . . . , Un^i) G N^~^ and (i G N, define 

= (mi,...,m„_i,c/) e N". 

Lemma 3. For any fi G N"~^ and (i G N, (/i, (i) G B{I) if and only if there exists 
k, where 1 < k < s, such that dk-i < d < dk and fi G B{Ik-i). 

Proof. {li,d) G B{I) if and only if there is no m G Min(J) such that m\X^^''^\ As 
dk-i d < dk, we only need to see that there is no m G Min(J) with deg^,^ m < dk_i. 
But this is equivalent to requiring that // G B{Ik^i). □ 

For a set of vectors U and an integer d, define 

U ®d = {{u,d) -.ueU}. 
Theorem 4. Irr(I) = IJfc=i (-^''"K-^fc-i) \ -^''"K-^fc)) ® '^k, which is a disjoint union. 

Proof. Assume fi G Irr(/fe_i) \ Irr(/fc). We first show that {fi,dk) 1 G B{I) and 
fiQl G i?(/fc„i)\S(Jfc). Since /i G Irr(4_i), we have /iBl G B{Ik-i), so {fi,dk)Ql = 
{fiQ 1, (ifc ~ 1) G B(/) by Lemma 3. Also, by Lemma 2, there is no 7 G -B(Jfc„i) such 
that 7 > yU0 1, in particular no 7 G such that 7 > /i0 1, as C 

Thus /i 1 ^ B{Ik), otherwise we would have /i G Irr(/fc) which contradicts the 
assumption on fi. 

For (/i, dk) G Irr(/), we need to prove that (/x, cifc) 1 is maximal in B{I). Assume 
otherwise, say (7,0?) G B{I) and (7,6?) > {fi,dk) 1. Then d > dk or d = dk — I. 
If d > dk, then 7 G B{Ij) where A; < j < s by Lemma 3. Since 7 > /U. 1 and 
B{Ik) is a 5-set, 7 G implies 1 G B{Ij) C too, a contradiction. If 

d = dk — I, then 7 > fi Q 1. Note that (7, dk — 1) ^ B{I) implies 7 G B{Ik-i) by 
Lemma 3. However, G Irr(Jfc_i) so there is no 7 G B{Ik-i) such that 7 > /i 1, a 
contradiction. Hence such (7,0?) does not exist. Consequently, {fi,dk) G Irr(J). 

Conversely, assume {fi,d) G Irr(/), we need to prove that there exist some 1 < 
k < s such that d = dk and fi G Irr(/fe_i) \ Irr(/fc). By Lemma 2, {fi,d) G Irr(/) 
implies 

{fi,d)ei e B{i), (3) 

and there is no (7, /) G B{I) such that 

{^,i)>{fi,d)ei. (4) 

By Lemma 3, (3) implies there exists k such that /i 1 G B{Ik~i), and 

4-1 < - 1 < 4- (5) 

By Lemma 3 again, (/i 1,4 ~ 1) ^ Then (4) and (5) imply that d = dk- 

(4) and (5) also imply that there is no 7 such that 7 G B{Ik~i) and 7 > /i 1, so 
G Irr(4_i). 
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It remains to prove /i ^ Irr(Jfc). Assume fi G Irr(/fc). Then /i 1 G B{Ik). 
By Lemma 3, (// l,dk) G B(/) and (/i l,dk) > (/i, (ifc) Q 1, contradicting to 
(/i, 4) e Irr(/). Thus /i G Irr(4_i) \ Irr(4). □ 

Theorem 4 gives us the following recursive algorithm for finding irreducible de- 
composition of monomial ideals. Suppose we are given / = (X"^ , • • • , X'^p) and fixed 
variable order xi < ■ ■ ■ < x„. We encode the set {ai, . . . , ap} as a tree T of height 
n. Our algorithm Irr(T) takes T as input and produce Irr(/) as output. That is, 
Irr(J) = Irr(T). 

Recursive Algorithm: Irr(T) 

Input: T, a tree encoding Min(/) 

Output: S, a set (or a tree) representing Irr(/) 

Step 1. Start at the root of T. If the height of T is 1, then T consists of a few leaves; 

let d be the largest label on these leaves and let S := {d}. 

Return S (and stop the algorithm). 
Step 2. Now assume T has height at least two. Set S* := { }. 

Step 3. Suppose do < di < ■ ■ ■ < dg are the labels of the children under the root of T, 
and let % be the subtree extending from dk, < k < s. 
Note that the root of is the node labeled by d^, but now unlabeled. 
Find Vo := Irr(7o) by recursive call of this algorithm. 
For k from 1 to s do 

3.1. Find % := MinMerge(7fc_i, 7^), and delete Tk-i. 

3.2. Find Vk := Irr(7fc) by recursive call of this algorithm. 

3.3. Find V := Vk-i \ Vk, delete 14-i, and S := Merge(^, V" 4)- 
Step 4. Return (S). 

Example 5. We end this section by demonstrating how the algorithm is used to 
decompose the ideal I = {x^,y^,x^y'^z'^,xy^z'^,x'^yz^). First represent the monomials 
as a tree with variable order x < y < z, where Tk 's are the subtrees extending from 
the node with label dk, /c = 0, 1, 2, 3. 

Figure 4-5 show the process of finding the irredundant irreducible decomposition of 
I. For each Tk, inductively MinMerge the subtrees from left to right, corresponding to 
Step 3.1 in the Recursive algorithm. See Figure 4- In Figure 5 we call the procedure 
Irr{ ) for each Tk to compute Irr{Tk), corresponding to Step 3.2. Since the height 
of Tk is 2, we bind each leaf that is not in the most-right side of Tk with the node 
of height 2 on the next path - just do the shifting in adjacent paths, see Figure 5. 
Finally we find the paths in Irr{Tk-i) that are not in Irr{Tk). The one with a mark 
X in Irr{Tk) is discarded. Then bind the resulting paths with dk. The irreducible 
components can be read from the last figure: 

Irr{I) = {(4, 4, 2), (4, 2, 3), (3, 3, 3), (4, 1, oo), (2, 3, oo), (1,4, oo)}. 
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Figure 5. Shifting step. 
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6. Incremental Algorithm 

In this section we shall present an incremental algorithm based on the idea of 
adding one generator at a time. This algorithm can be viewed as an improvement 
of Alexander Dual method ([14, 17]). We maintain an output list of irreducible 
components, and at each step we use a new generator to update the output list. In 
[17], it is not clear how to select good candidates that need to be updated, and the 
updating process there is also inefficient. Our algorithm avoids these two deficiencies. 
We establish some rules that help us to exclude many unnecessary comparisons. 

Monomial ideal are much simpler than general ideals. The next theorem tells 
us that monomial ideals satisfy distribution rules for the operations "+" and "fl" . 
These rules may not be true for general ideals. 

Theorem 6 (Distribution Rules). Let Ji, . . . , /f, J be any monomial ideals in K[X]. 
Then 

(a) + . . . + It) n J = hn J + . . . + Itf] J, and 

(b) (Ji n . . . n Ji) + J = (/i + J) n . . . n (/* + J). 

Proof. By induction, we just need to prove the case for t = 2. Note that (b) follows 
form (a), as 

(Ji + J) n (/2 + J) = /in(/2 + J) + Jn(/2 + J) 
= Ji n /2 + /i n J + J n /2 + J 
= hnh + J. 

To prove (a) for the case t = 2, suppose h is a generator for (Ji + I2) H J. Then 
h must be in (Ji + I2) and J. Since {Ii + I2) r\ J is also a monomial ideal, h is a 
monomial. The fact that h E Ii + I2 implies that h is in either Ii or l2- Hence h is 
in Ji n J or in /2 n J, so /i G /i n J + /2 n J. Going backward yields the proof for 
the other direction. □ 

Theorem 6 gives us an incremental algorithm for irreducible decomposition of 
monomial ideals. Precisely, we have the following situation at each incremental 
step: Given the irreducible decomposition Irr(J) of an arbitrary ideal / and a new 
monomial X" where a = (oi, . . . , a„) G N", we want to decompose 1 = 1+ {X"'). 
By the distribution rule (b), 

7=1 Pi m'^j+(X°)= fl (m^ + (X")). (6) 

\/3eIrr(/) / /3elrr(/) 

We need to see how to decompose each ideal on the right hand side of (6) and how 
to get rid of redundant components. We partition Irr(J) into two disjoint sets: 

Tf = {/5 e Irr(/) : a 7^ (3}, and (7) 
= {/? G Irr(/) : a ^ (8) 
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Note that if X" e / then = 0. For each /9 G Tf , we have X° e m^, thus 

+ (X°) = m^. (9) 
For each f] G T^, we have X" ^ m^. In this case, we spht (X°) as 

n 

(x-) = f|(x;o. 

By the distribution rule (b), we have 

71 

i=i 

Define 

/^("J) = . . . , bj^i, aj, bj+i, ...,bn), I <j <n. 
Since a -< /3, we have aj < bj for all 1 < j < n. Hence + (x"^) = m^^"'^\ and 

n 

+ (X") = Pi m^^'^''^ (10) 
i=i 

Therefore, 

Irr(J) = MaxMerge (Tf , {/J^"'^) : P e and 1 < j < n}) . (11) 

It remains to see which of the components in the right hand side of the above 
expression belong to Irr(/), so others are redundant. 

Lemma 7. Tf C /rr(7). 

Proof. Let Pi G . By equation (11) if Pi ^ Irr(J), then there exists some P2 G 
such that Pi is maxmergeed by /?2°''''' for some j, ie. Pi < p'^'^\ Since /52°''''''' < P2^ 
Pi < /52"'"'^ implies that Pi < P2, which contradicts with the fact that Pi,P2 G Irr(/). 
Hence Pi G Irr(/) as claimed. □ 

Lemma 7 shows that the elements in will be automatically in Irr(/). Now we 
turn to the components P^°''^\ For P & T2, define 

Mf3 = {m G Min(/) : m\X^}. (12) 

For m G Afg, if deg^,^ m = bu, then we say m matches p in It is possi- 
ble that one monomial matches P in multiple variables. For example, with / = 
{x^,y^, z'^,xy,xz,yz) and P = (1,1,2) G Irr(J), the monomial xy matches P in x 
and y. We say m matches P only in Xu if deg^,^ m = bu and deg^,^ m < bk for all 
k ^ u. 

Lemma 8. For each P = (61, . . . , 6„) G Tl* and each 1 < u < n, there exists m G Afg 
such that m matches P only in x„. 
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Proof. Note that a vector 7 G B(/) is maximal if and only if X'^ ■ x„ G / for every u. 
Since (3 G Irr(/), /50l is maximal in B{I). Thus, for each 1 < u < n, X^'^^-Xu G /, so 
there exists a monomial say m G Min(J) such that m|X^®^ ■ Then deg^^ m < bk 
for k ^ u. If deg^^ m < 6„ as well, then m|X^®^, which implies that X^®^ G /, a 
contradiction. Therefore deg^^ m = bu- Note that X^®-^ ■ Xu\X^, so m G Afg. □ 

For any set of monomials A C IK[X], define max(A) be the exponent 7 such that 
XT = Lcm(A). 

Lemma 9. max(M^) = j3. 

Proof. By the definition of M^, we know that max(M^) < f5. By Lemma 8 we have 
max(M^) > p. Thus max(M^) = p. □ 

For k ^ u, let 

M, /c) = min{deg^^ m : m E Mp matching (3 only in x^}. (13) 
Note that d{[3, u, k) < bu- Define 

d{(3,u) = ^^max {d{p,u,k)}. 

Lemma 10. For each (3 G and 1 < u < n, g Irr{I) if and only if 

d{(3,u) < ttu. 

Proof. Suppose d{P,u) < a„. We want to prove that /J^"'"") G Irr(/). By Lemma 
2, this is equivalent to proving that p^°'''^^ 1 G B(J) and is maximal. Assume 
p{a,u) Q I ^ ^^jy rjnj^g^ ^^^^^ g^jg^g ^ ^ Min(/) U {X"} such that m|X^'"'"'ei. 

First note that m 7^ X"" because X" can not divide Thus m G Min(J), 

which implies X^*"'"''®^ G /. Since /5*^"'")0l < /50l, we have X^®-^ G /, contradicting 
to p e Irr(J). Jience 1 G B(J).^We next need to prove that 1 is 

maximal in B(/), that is, X^^"'"''®^-Xfc G / for every k. In the case for k = u,we have 
For any k ^ u, let m be any monomial in (13) such that deg^,^ m = 
d{P,u,k). Then deg^^ m = d{P,u,k) < d{P,u) < Ou, hence m|X^'"'"''®^ ■ Xk as 
degj,^ m = bk and deg^^. m < bj — 1 for j ^ u, k. 

Conversely, suppose g Irr(J). We want to prove that d{P,u) < Ou- We 

know that 1 is maximal in B(/). Thus X^^"'"'®-^ ■ Xk E I for every k. For 

any /c 7^ u, suppose ' ®^ ■ Xfc is divisible by m G Min(J) U {X"}. Then 

deg^^ m<au-l<bu, deg^, m<bj-l, j ^u,k, (14) 

and deg^^m < bk- As X^'"'"'©i G B(7) C B(/), m can not divide x^'"'"'®^. Hence 
deg^,^ m < bk- So m matches P only in Xk- Note that m 7^ X°, so m G M and thus 
m G Mg. It follows that d{P,u,k) < a^i — 1 by (14). Therefore, d{P,u) < as 
desired. □ 
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By the above lemma, for each (3 G we only need to find Mp and d{(3,u), 
which will tell us whether G Irr(J). This gives us the following incremental 

algorithm. 

Incremental algorithm 

Input: M, a set of monomials in n variables 

Output: Irr(/), the irredundant irreducible components of the ideal / generated by M. 
Step 1. Compute MinMerge(M) and sort it into the form: 

MinMerge(M) = {xl\ . . . . . . 

where Cj can be 00 and {X"^, . . . , X"'''} are sorted in lex order with variable 
order xi < . . . < x„. Set 

T := {(ai, . . . ,a„)}. 
Step 2. For each k from 1 to p do: 

2.1. Set the temporal variables V = ^ and a := a^- 

2.2. For every j3 E T with a 13 do 

V ■.= VU{(3]. 

2.3. For every (3 eT with a -< (3 do, 

• find Mp as defined in (12); 

• for 1 < M < n, compute d{[3, u), and if d{l3, u) < Qu then update 

V := 

2.4. Set T := V. 
Step 3. Output T. 

We next prove that there is a nice property of the above algorithm for generic 
monomial ideals, that is, the size of T is always non-decreasing at each stage when 
a new generator is added. This will allow us to bound the running time of the 
algorithm in term of input and output sizes. 

Theorem 11. Suppose I is generic and Min{I) = . . . , x^", . . . , 

where X"*-' 's are sorted in lex order with variable order Xi < . . . < x„. Let I = 
{xl\...,x^", . . Then \Irr{T)\ < \Irr{I)\. 

Proof. Keep notations as above. For every /5 G T^*, bn = Cn- Thus x^" is the only 
monomial in Mp that has degree in x„ larger than a^. Hence d{f3,n) < an and 
^(a.") £ Irr(/). By the equation (11) and Lemma 7, 

|Irr(/)| > |Tf I + : [3 G T^}\ = |Tf | + \T^\ = |Irr(J)|. □ 

The reader might wonder whether a similar statement holds in non-generic case 
as well. The answer is negative. Let / = {x^ , , z'^ , w"^ , x'^yz, xy'^w) C lK.[x,y, z,w] 
with lex order and x < y < z < w. Then 

Irr(/) = {(3, 3, 1, 1), (2, 3, 2, 1), (3, 2, 1, 2), (3, 1, 2, 2), (2, 2, 2, 2), (1, 3, 2, 2)}. 
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By adding X" = xyzw, we can see {3 = (2,2,2,2) e T^. Note that = 
{x'^yz,xy'^w, z'^, w'^}. Since d{P,u) = 1 = a„ for m = 1,2,3,4, no new P^""'^^ will be 
generated. Thus the number of irreducible components decreases by 1 instead. 

We find the irreducible components for the monomial ideal in Example 5 again 
by the flow of our incremental algorithm. 

Example 12. Decompose 

I = {x\ y\ xYz^ xyh\ x^yz^). 

Note: "/ " means 6 Irr(I) for corresponding /5, a and u, while "x " means 

not. 

Step 1. M = {x^,?/^,^~,xV^^a;?/=^^^xV^}■ 5ei T := {(4, 4, oo)}. 
Step 2. (i) For a = (3, 2, 2) do: 

2.1. V := <p. 

2.2. Since a -< (4,4, oo), V := 0. 

2.3. Let (3 = {4,4, oo). We find = {x\y^} . 

So we have d{(3,l} = Q(/), 4/9,2} = Q(/) and 4/3,3} = 0^/;. 
Then V := {(3, 4, oo), (4, 2, oo), (4, 4, 2)}. 

2.4. LetT := V. 

(11) Fora = (1,3,2) do: 

2.1. V ■= 0. 

2.2. Update V by V := {(4, 4, 2), (4, 2, 00)}. 

2.3. a -< (3,4,00). 

Letp= (3,4,cx)). We find Mf} = {y\x^yh^}. 

So d{f3, 1} = 0^/ 4/3, 2} = 2(/) and d{(], 3} = 2^x 

Then V := {(4, 4, 2), (4, 2, 00), (1, 4, cx)), (3, 3, 00)}. 

2.4. LetT := V. 
(Hi) For a = (2, 1, 3) do: 

2.1. V := 4>. 

2.2. \/:={(4,4,2),(l,4,oo)}. 

2.3. a -< (4,2,00), and a -< (3,3, 00). 

• Letf3= (4, 2, 00) . We find Mp = {x^ x^y^z^} . 

So d{f3,l} = 3(x), 4/3,2} = 0^/) and d{(3,3} = 2(/ ). 
Then V := {(4, 4, 2), (1, 4, 00), (4, 1, 00), (4, 2, 3)}. 

• Letp= (3,3,00). Then Mp = {x^y^z'^,xy^z^}. 
4/3,1} = l(^), 4/3,2} = 2(x), 4/3,3} = 2(/). 

So V := {(4, 4, 2), (1, 4, 00), (4, 1, 00), (4, 2, 3), (2, 3, 00), (3, 3, 3)}. 

2.4. LetT := V. 

Step 3. Output T= {(4, 4, 2), (1, 4, 00), (4, 1, 00), (4, 2, 3), (2, 3, 00), (3, 3, 3)} 
= {(4, 4, 2), (4, 2, 3), (3, 3, 3), (4, 1, 00), (2, 3, 00), (1, 4, 00)}. 
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Some preprocess can be taken right before Step 2 to improve the efficiency of the 
incremental algorithm. For each n G {1, . . . , n}, we partition M into disjoint subsets 
such that the monomials in each subset have the same degree in x^. We then store 
these information, which requires memory complexity 0{n ■ p). For each /3 G Tg", 
we can find Mg by only checking the monomials in the subset with degree in 
variable x„ for every u. Note that for generic monomial ideals each subset contains 
a unique monomial. In this case Afg contains n monomials, and it can be found by 
0(77,) operations, instead of 0{p) operations by scanning through the whole input 
monomial set. 

7. Time Complexity and Conclusion 

We estimate the running time of our algorithms by counting the number of mono- 
mial operations (ie. comparisons and divisibility) used. Our recursive algorithm 
depends heavily on the number of distinct degrees of each variable. Let Sj be the 
number of distinct degrees of xj where j = l,...,n. Then the total number of 
merge of subtrees used by the algorithm is at most 11^=1 ^j- Since each subtree has 
at most p leaves(ie. p generators), each merge takes 0{p^) monomial operations. 
Hence the algorithm uses 0{p^ ■ 11^=1 "^j) monomial operations. This algorithm is 
more efficient for highly non-generic monomial ideals. The benchmark analysis in 
[20] compare several algorithms based on various slicing strategies, including our 
recursive algorithm. It is shown there that our algorithm performs as a very close 
second best one. 

The running time of our incremental algorithm is harder to estimate for general 
ideals. For generic ideals, however, we can bound the time in terms of input and 
output sizes. More precisely, suppose 

/ = (xr,...,<",x"s...,x°'') 

is a generic monomial ideal in K[X] where X^'^'s are sorted in lex order with variable 
order xi < . . . < x„. For < A; < let 

/(fc) = (x?,...,CX°S...,X"^). 

All these ideals are generic. By Theorem 11, we have 

1 = |Irr(/(o))| < |Irr(/(i))| < ■•■ < |Irr(/(,))| = |Irr(J)|. 

In an arbitrary stage of the incremental algorithm, we try to find the irreducible 
components of /(fc) from those of I(k-i)- For each j3 G Irr(/(fc_i)), only those P in T^*" 
(as defined in (8)) need to be updated. Note that I is generic, by the preprocess 
Mj3 can be found in 0{n) operations. The numbers d{f3,u, k), 1 < u ^ k < n, can 
be computed by scanning through the monomials in once, thus using only 0(n) 
monomial operations. Then the numbers d{(3,u), 1 < u < n, can be computed in 
O(n^) operations. Hence for each /3 G T2'', Step 2.3 uses at most 0{n + n?) = O^n^) 
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monomial operations. Since T D T^*" has at most i elements where i = |Irr(J)|, 
Step 2.3 needs at most 0{n'^£) monomial operations. Therefore, the total number 
of monomial operations is at most 0{n'^pt). In fact, T^* is usually a small subset of 
T, the actual running time is much better than our worst-case estimate indicates. 

We also want to point out that for generic monomial ideals, the incremental 
algorithm is an improved version of the recursive algorithm. Suppose we add the 
new monomial X°''' into I(k~i)- In Step 3.2 of the recursive algorithm, we need to 
compute Irr(7^). But in Step 2.3 of the incremental algorithm, only /3 G T^* need 
to be updated. We have the observation that is a small subset of Irr(7fc) ® c„. 
By this observation we conclude the incremental algorithm is more efficient than the 
recursive algorithm for generic monomial ideals. In non-generic case, the comparison 
is not clear. 

In all previous algorithms (including our recursive one) for monomial decompo- 
sition, the storage in the intermediate stages may grow exponentially larger than 
the output size. Our incremental algorithm seems to be the first algorithm for 
monomial decomposition that the intermediate storage is bounded by the final out- 
put size. Note that the output size i can be exponentially large in n. In fact, it is 
proven m [1] that^ = 0(p[tl) for large p. Since the output size can be exponential in 
n, it is impossible to have a polynomial time algorithm for monomial decomposition. 
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